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The spin-Hall effect in the two-dimensional electron gas (2DEG) generates symmetric out-of-plane spin Sz ac- 
cumulation about the cuiTent axis, in the absence of external magnetic field. Here we employ the real space 
Landauer-Keldysh formalisml l) by considering a four-terminal setup to investigate the circumstances in which 
this symmetry is broken. For absence of Dresselhaus interaction, starting from the applied out-of-plane B corre- 
sponding to Zeeman splitting energy A = to 0.5 times the Rashba hopping energy t^Q, the breaking process 
is clearly seen. The influence of the Rashba interaction on the magnetization of the 2DEG is studied herein. For 
coexisted Rashba t^Q and Dresselhaus t^Q spin-orbit couplings in the absence of B, interchanging t^Q and 
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reverses the entire accumulation pattern. 



PACS numbers: 73.21 .-b, 71.70.Ej, 75.60. Ej 
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The generation and transport of spin currents dominates the 
appHcations of spintronics. The spin-orbit (SO) interaction, 
which couples the electric degree of freedom with the mag- 
netic one serves as the mechanism to achieve this. A number 
of basic designs' ^* for spintronic devices, such as field-effect 
switches, spin transistors, spin filters, and spin waveguides, 
have been proposed by taking advantage of this interaction 
to control spins. One of the phenomena originating from the 
SO interaction is the spin-Hall (SH) effect' in which a trans- 
verse spin current is induced by 3. longitudinal electrical cur- 
rent. The semiclassical SO force*^ proportional to (p x e^)^, 
oppositely deflects the spin up (cr = 1) and down (a — —1) 
wave packets with momentum p in the transverse directions 
so that different spins accumulate in the lateral edges. The 
effect is particularly notable in that it produces a spin current 
with no magnetic field applied, and no accompanying charge 
current present, rj^pp^t experimental work using scanning 
Kerr microscopy'^'^'^'g in n-type unstrained GaAs, strained 
InGaAs'^ and 2DEG,'^ have inspired a host of theoretical 
studies on SH effect. 

The 2DEG confined in the InGaAs/InAlAs semiconductor 
heterostructure''" possesses intrinsic SO coupling known as 
the Rashba and Dresselhaus interactions. The inversion asym- 
metry of the structure gives this system the Rashha-interac- 
tion with adjustable strength via the gate voltage"^^ while 
the bulk inversion asymmetry' gives rise to the Dresselhaus 
interaction whose strength is material dependent. A num- 
ber of studies regarding the intrinsic SH (ISH) effect in the 
2DEG have been reported, but most of these focus on evaluat- 
ing the spin current which is, due to its nonconservation, not 
easily measured. The SO coupling leads an electron to pre- 
cess around a momentum-dependent effective magnetic field 
which produces a source or a sink of the spin current in the 
continuity equation. Due to tWsMionconservation the spin cur- 
rent is not uniquely defined. ^^l^ Although a possible defini- 
tion of a conserved spin current has been suggested in Ref. [3, 



the present work follows Ref. T by employing the Landauer- 
Keldysh formalism in the four-terminal setup to investigate a 
more directly measurable physical quantity: the accumulation 
of the out-of-plane spin Sz- If only the Rashba interaction ex- 
ists then the accumulation induced by the electric potential de- 
posits spins symmetrically along the transverse direction, i.e., 
the up-spins accumulate on one side while the down-spins ac- 
cumulate symmetrically on the other. The majority of current 
studies focuses on how this SH symmetry (SHS) is generated. 
By contrast, in the paper we answer the essential question of 
the circumstances under which the SHS is brokep- 

With the help of the finite difference method,''^ the linear 
Rashba and Dresselhaus model with applied field B is, in the 
tight binding limit, expressed as 
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where cl^^. (Cm'a') denotes the creation (annihilation) op- 
erator for a = +1 (spin-up) or cr = —1 (spin-down) at 
the site m = {m.x,my). The electric potential and disor- 
der can be accounted for by the on-site energy Em- To see 
how SHS is destroyed either by B or coexisted Dresselhaus 
and Rashba interactions, we consider the clean limit and set 
the tight-binding bottom energy, for convenience, to be — 4to 
(corresponding to Em = 0), with to = I2ma^ being the 
hopping energy. The Zeeman sphtting A = —ehB/mc is 
induced by the magnetic field while the Rashba (Dressel- 

a{(3)/2a is taken into account 



haus) SO coupling t^Q^"* 



by the neareast-neighbor hopping matrix element t 
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the identity matrix in spin space is Ig, the lattice constant is 
a, and the Rashba (Dresselhaus) coupling strength is a (/?). 
In order to subject the conductor to an electric potential dif- 
ference along the x axis we consider the Landauer setup with 
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the four ideal leads p = 1 (left), 2 (right), 3 (bottom), 4 (top) 
shown in Fig. [Tfa). 

To acquaint the reader with the method employed herein we 
present here a brief review of the Landauer-Keldysh Formal- 
ism. In general, in a conductor the non-equilibrium spin accu- 
mulation at time t, {S^{t)) = fi/2j2^a{cl,^{t)cma{t)) = 
;i/2X;„(T[(mcr|G<('=)(t,t)|mcr)?i/i] , depends on the 
switching time ts (at which leads are brought into contact) 
via the lesser Green function 

G<'\h,t[) (2) 

/OC pOO 

which includes both the steady state (second line), and also the 
transient state (third line). For the measuring time t much later 
than tg (this is the case of interest here) one can approximately 
write in Eq. ^ tg = — oo. The transient solution can then be 
neglected, and G<(^^(fi, = G^''^^(t) depends only on the 
the time interval t = ti — t[. A Fourier transformation of Eq. 
then yields the kinetic equation 

G<'-''\E) = G'-(=)(i;)S<('=)(£;)G'-('=)t(^) (3) 
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p = 4 



the retarded Green function G'"(^^(i?) = 
1 thot fijg steady 
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with 

[E - H^'=^ - so that 

accumulation is expressed as S^{t 

The lead interacts with the conductor through the self- 
energyil^) s'^^ — J2p'^i^\ ^^^^ matrix elements 

(m,a|s(^V,^') ^ s'lX' = 4g<^^Hrm.r'^,)S^^, 
for m and m' (in the conductor) being adjacent points to rm 

t(c)(T(t' 



-pm m' *-* Otherwise. The hopping 
energy for these points allows electrons to flow through the 



and r'^i (in leads), and 
energy for these points £ 
interfaces. While the lesser self-energy is written as 

S<(=)(£;) = -2ilm'Ep{E - eVp)f{E ~ eVp), 



(4) 



where f{E) = 1/ [1 + exp(-£'/fcBT)] is the Fermi-Dirac 
distribution and the quasi-particle escaping time, in the con- 
ductor, is inversely proportional to Im S*^^^ (E). The retarded 
Green function g^P^ for the isolated lead p can be evaluated by 
considering the single particle Green operator g^^*^ — 1/{E + 
iO+ — H'-p^) in the eigenfunction expansion g'^P') {ri,r2) = 
'En^^^rF^{rii,ru)^i''^*{r2i,r2t)/{E + iO+ - el^h where 
r = (n, r() is the position vector (within the lead) composed 
by the longitudinal component r; and the transverse compo- 
nent rt with n and 7 accounting the transverse and longitu- 
dinal modes, respectively. The eigenfunction ip^j {ri, rt) ~ 
2 sin(n7rrf /W) sin(77rn / L)/ {a\/WL), with the normaliza- 



tion Yln-^ V-'n7 (r')V'ri7 (r) = 5rr' , IS obtained by solving 
the eigenequation Hpipn-^ {ri,rt) = [—{h^ /2m)d'^ / drf + 
ycoif{Tt)]iJn^ {ruTt) = en'y'ipnjiri,rt) under the hard-wall 
boundary condition in which the confining potential 
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FIG. 1 (Color online) (a) Illustration of the Landauer four-terminal 
setup with the out-of-plane B. Contacting the 2DEG (green square) 
with four ideal semi-infinite leads (yellow rectangles) p = 1 (left), 
2 (right), 3 (bottom), 4 (top) which are biased with eVi — —eV2 = 
eV/2, eVa = eVi = and eV = lO^^to, the spatial out-of-plane 
5^ accumulations, in units of ft/2 with Zeeman splitting A being 
varied from to O.StgQ, are plotted in (b). The magnetization, ob- 
tained by summing over every accumulation on each site, in units of 
h/2, decreases with t^Q but increases with A as shown in (c). 



00 (0) outside (inside) the lead p. The width (length) of the 
lead is W (L). For semi-infinite lead L —> 00, g(^')(ri,r2) 
can be directly computed by replacing the summation of lon- 
gitudinal modes with d-jL/TT. 

Consider now the InGa As/In Al As heterostructure''^' with 
typical parameters. The effective mass m = 0.05me (nie is 
the electron mass) and the lattice constant a — 3 nm yield 
the hopping energy to = 84.68 meV. We set tgp = O.Uq, 
eVi = -eV2 = eV/2, eVs = eVi ^ 0, eV = 10-'% 
and conductor size to be 8a x 8a. Select the Fermi energy 
Ep — — 3.8t close to the band bottom at — 4io, so that the 
tight-binding approximation valid. To examine how A and the 
coexistence of t^Q and tgp affect the SHS, we show the spa- 
tial 5"^ accumulation in units of h/2 in Fig. [T](with t^Q — 0) 
and Fig. in (with A = 0). 

Obviously, the applied field B polarizes the 2DEG or, from 
perspective of the band theory, it induces a Zeeman splitting 
A such that the SHS is destroyed. Starting from A = and 
going to A = 0.5tgQ, two effects are found to break the sym- 
metry [see Fig. [Tfb)]: (i) The area of the majority spins is 
enlarged, (ii) The magnitude of the majority magnetization 
is increased. On the other hand, for parameters A and t^Q 
varying between and O-Uq, the magnetization (or the total 
z-polarization) of the system, obtained by summing over ev- 
ery Sz on each site, is decreased with increasing Rashba inter- 
action as shown in [Fig. [Uc)]. To explain this, we note that. 
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FIG. 2 (Color online) (a) The spatial Sz accumulation with 
{tsoy^so) ~ (O.lto, O.OSto) ■ (b) The Sz accumulations in units 
of h/2 as functions of positions with parameters {tsoy^so) ~ 
(O.Uo,0.05to) (green) and {tsc^so) = (0.05to, O.Uo) (blue). 
The red line sums up these two functions. 



an accumulation pattern different from the one which results 
from pure Rashba interactions. In particular, destructions of 
the SHS are found in two special cases: in the presence of a 
magnetic field B and in the presence of coexisting Rashba 
and Dresselhaus couplings. In the former case, beginning 
with both A = and t^p = 0, the applied field B breaks 
the SHS by not only extending the area of the majority spin 
accumulation, but also by strengthening its magnitude. Mean- 
while, the gate-voltage-tunable Rashba interaction reduces the 
magnetization induced by the Zeeman splitting. In the latter 
case (where both SO couplings are present), interchanging the 
Rashba and Dresselhaus interactions reverses the whole accu- 
mulation pattern. These features thus provide an electric con- 
trol of a magnetic property. 

This work is supported by the Republic of China National 
Science Council Grant No. 95-2112-M-002-044-MY3. 



under the Rashba interaction, electrons precess with respect 
to an in-plane (x-y) effective Rashba magnetic field which 
yields vanishing mean z-polarization so that magnetization is 
reduced. 

Even B is absent, the SHS can also be broken by the co- 
existence of tg^Q and tgQ. Consider the Dresselhaus coupling 
tfo = O.lio (0.05io), corresponding to /3 w j {kj) = 0.5 
(0.025) eVA, assuming a z direction quantum pwell width 
22 A (31 A) with typical value of the coefficient"® 7 25 
eV A . Our results indicate that the SHS is preserved by the 
presence of tgp or tgp alone, but it is destroyed if the Rashba 
and Dresselhaus interactions coexist. Figure |2f a) shows the 
accumulation pattern for t^Q — Q-Uq, tgQ ~ O.OSto, and 
S = 0. In comparison to the A = result of Fig. [Hb) it is 
tilted upwards in the left hand region, and downwards in the 
right hand region. Furthermore, if the strength of the Rashba 
tgQ and the Dresselhaus t^p interactions is interchanged then 
the pattern is entirely reversed, i.e., the up (down) accumula- 
tion becomes the down (up) accumulation. To illustrate this 
effect we label the spin component 5*^ by numbering the po- 
sition of each site row by row. For example, m — (a, a) is 
labeled by position 1 and m = (8a, 8a) is labeled by position 
64. The accumulation of the component Sz in the absence of 
the field B is plotted in Fig. |2jb) as a function of position for 
two special cases: (i) {t§o^tso) = (0.1to,0.05<o) denoted 
by the green line, and (ii) [tso^^so] ^ (0.05to, O.lto) de- 
noted by the blue line. Summing up these two functions at ev- 
ery position we obtain zero accumulation everywhere (the red 
line). This suggests an inverse correlation between the cases 
(i) and (ii). In the case of equal strengths, t'gQ — tgQ, one 
therefore expects that Sz accumulates nowhere, since swap- 
ping tg^Q and tgQ has no effect. Finally, we recall that we 
address a finite 2DEG. Comparing our results with previ- 
ous works on infinite systems we identify the predicted sign 
change' in the spin-Hall conductivity. 

In conclusion, we investigate the out-of-plane Sz spin ac- 
cumulation using the Landauer-Keldysh formalism in a four 
terminal setup. Taking into account the Rashba tgQ and the 
Dresselhaus t^Q couplings and Zeeman splitting A we obtain 



References 

[1] B. K. Nikolic, S. Souma, L. P. Zarbo, and J. Sinova, Phys. Rev. 

Lett. 95, 046601 (2005); B. K. Nikolic, L. R Zarbo, and S. 

Souma, Phys. Rev. B 73, 075303 (2006). 
[2] S. Datta and B. Das, Appl. Phys. Lett. 56, 665 (1990); J. Nitta 

et ai, Appl. Phys. Lett. 75, 695 (1999); Francisco Mireles and 

George Kirczenow, Phys. Rev. B 64, 024426 (2001); T. Koga et 

al., Phys. Rev. Lett. 88, 126601 (2002); X. F Wang et al, Phys. 

Rev. B 65, 165217 (2002); D. Frustaglia and K. Richter, Phys. 

Rev. B 69, 235310 (2004); J. Carlos Egues, Guido Burkard, D. 

S. Saraga, John Schliemann, and Daniel Loss, Phys. Rev. B 72, 

235326 (2005). 

[3] J. Sinova, D. Culcer, Q. Niu, N.A. Sinitsyn, T. Jungwirth, and 
A.H. MacDonald, Phys. Rev. Lett. 92, 126603 (2004). 

[4] Branislav K. Nikolic, Liviu P. Zarbo, and Sven Welack, Phys. 
Rev. B 72, 075335 (2005). 

[5] Y.K. Kato, R.C. Myers, A.C. Gossard, D.D. Awschalom, Sci- 
ence 306 1910 (2004). 

[6] J. Wunderlich, B. Kaestner, J. Sinova, and T. Jungwirth, Phys. 
Rev. Lett. 94, 047204 (2005). 

[7] V. Sih, R.C. Myers, Y.K. Kato, W.H. Lau, A.C. Gossard, D.D. 
Awschalom, Nature Phys. 1 31 (2005). 

[8] N. P Stem, S. Ghosh, 1 G. Xiang, M. Zhu, N. Samarth, and D. 
D. Awschalom, Phys. Rev. Lett. 97, 126603 (2006). 

[9] J. Nitta, T. Akazaki, H. Takayanagi, and T. Enoki, Phys. Rev. 
Lett. 78, 1335 (1997). 
[10] B. Das, D.C. Miller, S. Datta, R. Reifenberger, W.P Hong, PK. 
Bhattacharya, J. Singh, and M. Jaffe, Phys. Rev. B 39, R1411 
(1989). 

[11] G. Dresselhaus, Phys. Rev. 100, 580 ( 1955). 

[12] P Zhang, J. Shi, D. Xiao, and Q. Niu,'cond-mat/050350 5J 
[13] S. Murakami, N. Nagaosa and S. C. Zhang, cond-mat/03 1 0005! 

to appear in Phys. Rev. B. 
[14] J. Shi, P Zhang, D. Xiao, and Q. Niu, Phys. Rev. Lett. 96, 

076604 (2006). 

[15] Supriyo Datta, Electronic Transport in Mesoscopic Systems 

(Cambridge University Press, 1997). 
[16] W. Knap, C. Skierbiszewski, A. Zduniak, E. Litwin- 

Staszewska, D. Bertho, F. Kobbi, and J. L. Robert, Phys. Rev. 

B 53, 3912 (1996). 
[17] N. A. Sinitsyn, E. M. Hankiewicz, W. Teizer, and J. Sinova, 

Phys. Rev. B 70, 081312(R) (2004). 



